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We construct gravitating superconducting string solutions of the i/(l)iocai x f (l)giobai model solv- 
ing the coupled system of Einstein and matter field equations numerically. We study the properties 
of these solutions in dependence on the ratio between the symmetry breaking scale and the Planck 
mass. Using the macroscopic stability conditions formulated by Carter, we observe that the cou- 
pling to gravity allows for a new stable region that is not present in the flat space-time limit. We 
match the asymptotic metric to the Kasner metric and show that the relations between the Kasner 
coefficients and the energy per unit length and tension suggested previously are well fulfilled for 
symmetry breaking scale much smaller than the Planck mass. We also study the solutions to the 
geodesic equation in this space-time. While geodesies in the exterior space-time of standard cosmic 
strings are just straight lines, test particles experience a force in a general Kasner space-time and 
as such bound orbits are possible. 

^ | PACS numbers: 11.27.+d, 98.80.Cq, 04.40.Nr 
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I. INTRODUCTION 



Particle physics theories beyond the Standard Model generally predict phase transitions in the early universe, during 
0-| which topological defects can appear [l[. Since they are, by definition, topologically stable-they may survive up to now 
and have detectable effects. In particular, one-dimensional defects, called cosmic strings @, appear quite generally 
and are thought to reach a scaling regime, so that their contribution to the energy density of the universe remains finite 
and constant. Due to the fact that these objects can be extremely heavy they were believed to be a possible source of 
the density perturbations that led to structure formation and the anisotropies in the cosmic microwave background 
(CMB) However, the detailed measurements of the CMB power spectrum as obtained by COBE, BOOMERanG 
and WMAP demonstrated that cosmic strings cannot be the main source for these anisotropies. However, in recent 
years it has been suggested that cosmic strings should generically form at the end of inflation in inflationary models 
resulting from String Theory Q such as brane inflation Q. Moreover, cosmic strings seem to be a generic prediction 
of supersymmetric hybrid inflation [6| and grand unified based inflationary models [7j. Even though the origin of 
these cosmic superstrings is String theory, their properties can be investigated in the framework of field theoretical 

models [Hill]. 

_^ . In general cosmic strings can not end. This allows for two kinds of strings: loops and infinitely extended strings. 
Infinite strings arc thought to be relatively straight on macroscopic scales and their width is typically much smaller 
than their extent. A string can thus be described in good approximation by a one-dimensional object, characterized by 
quantities integrated over a plane orthogonal to it. This procedure is generally not well-defined in curved space-time, 
but can be used here since the metric generated by the string is "nearly Minkowskian" reasonably far away from it. 
For infinite straight strings, there are in fact only two relevant quantities: the energy per unit length U and the tension 
T. In the simplest field theoretical model, namely the Abelian-Higgs model, these quantities are equal This 
can be related to the necessary Lorentz invariance along the string axis. But if neutral currents or some microscopic 
structures (for instance wiggles) are taken into account [13| . then the energy per unit length can be larger than the 
tension. If the tension would be larger than the energy per unit length the string would be unstable under transverse 
perturbations. The relation between the energy per unit length and the tension of a superconducting string solution 
of the ?7(l)i oca i x J7(l) g iobal model has been discussed in detail in [2 using the formalism developed by Carter [l5l.ll6j 
and it has been suggested that the equation of state is of logarithmic form. This has been confirmed numerically in 
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It is also of interest for a possible observation of cosmic strings how the space-time of such an object looks like and 
how test particles would move in this space-time. In the case of infinitely thin standard cosmic strings, i.e. cosmic 
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strings without additional structure that fulfill U = T, the space-time is locally flat and geodesies are just straight 
lines. However, the space-time is globally conical with deficit angle A = 8irGU and as such light gets deflected by 
a cosmic string. Next to the observation of cosmic string signals in the power and polarization spectra of the CMB 
p^| - [20j this lensing property of cosmic strings has been suggested to be the prime signature of these objects. Taking 
the finite width of the cosmic string into account [2l[ bound orbits of massive test particles are possible, while massless 
particles can only move on escape orbits. 

The exterior metric of a superconducting string carrying timelike and spacelike currents, respectively, has been first 
discussed in [22J, while lensing properties of a cosmic string with a lightlike current have been investigated in 23 1. 
The space-time of cosmic strings with a non-degenerate energy-momentum tensor has been given in [24| and it has 
been shown that the parameters in the general Kasner space-time [25j can be given in terms of U and T. This will 
describe e.g. the exterior of superconducting strings as well as strings with wiggles [2a |. 

When considering concrete field theoretical models of gravitating cosmic strings to describe their microscopic 
properties, we need to solve the set of coupled Einstein and matter field equations numerically. This has been done 
for Abelian-Higgs strings without currents in [UHl]. In this paper, we are interested in solving the full set of coupled 
Einstein and matter fields equations of a E/(l)i OC ai x C/(l)global model describing superconducting string solutions with 
either timelike or spacelike currents in curved space-time. We will hence be able to determine the metric functions 
on the full interval from the string axis out to infinity. 

Our paper is organised as follows: in Section II, we discuss the field theoretical model describing gravitating 
C^(l)iocai x t^(l)giobai superconducting strings. In Section III, we discuss our numerical results and in particular match 
our solutions to the Kasner solutions. In Section IV, we comment on test particle motion in these space-times and 
we conclude in Section V. 



II. THE MODEL 



In the following, we will consider the £7(l)i OC ai x ^(l)giobal m °del in curved space-time. The action reads 

S= [ d i x^(-J—R + £ m ) (1) 



where R is the Ricci scalar and G denotes Newton's constant. The matter Lagrangian C m reads: 

c m = \d^(d^)* - -F^ + ^WO* - V(<f>, (2) 

with the covariant derivative D^<j> = V M — ieA^cj) and the field strength tensor = V^-A^ — V vA^ = d^A v — d^A^, 
of the U(l) gauge potential A^ with coupling constant e. The fields cf> and £ are complex scalar fields with potential 

v(<t>, = y (00* - vif + y£C (#* - H) + • (3) 

In cylindrical coordinates {t, r, 9, z} we choose the following Ansatz for the matter fields 

A^dx^ = - (n — P(r)) d9 , 0(r, 6) = r/ 1 h(r) exp(in^) , £(r, t, z) — r)if(r) exp(ifcz — itut) (4) 

and 

ds 2 = N 2 {r)dt 2 - dr 2 - L 2 (r)d9 2 - K 2 (r)dz 2 (5) 

for the metric. In the flat space-time limit, i.e. for G = we can define a Lorentz-invariant quantity w := k 2 — uj 2 
which is typically used to categorize the solutions into the timelike, spacelike and chiral type for w < 0, w > and 
w = 0, respectively. In the following, we will adopt the viewpoint that we can always go to a suitable frame of 
reference to choose either tu 2 = or k 2 = 0. Hence w — k 2 and w = — oj 2 corresponds to the spacelike and the timelike 
case, respectively. 

A. Equations of motion and boundary conditions 

The dynamics of the metric is given by the Einstein equations which read 

R pv = -8ttG (t^ - \g„ v T\ , (6) 
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where T = T£ is the trace of the energy-momentum tensor given by 

TZ = 6ZL m -2g^-^ . (7) 

The components of the Ricci tensor then are [13] 

JLKNJ R r__N^_L^_K^_ „e_ (NKL'Y _ _ (NLKj 

Kt ~ LKN ' Kr ~ N L K 1 Ke ~ NKL ' "* ~ NLK ' ( *> 

where here and in the following the prime denotes the derivative with respect to r. 
We use the rescalings 



r-> , L -> , (u,k) ->■ V^im(w.fc) , (9) 



where Mh = y/XiVi corresponds to the mass of the Higgs field. The components of the energy-momentum tensor and 
the field equations will then depend only on the following dimensionless coupling constants 

e 2 A 
a:=— , K :=87rG?7 2 , 7i :=-!,* = 2,3 . (10) 
Ai Ai 

The constant k corresponds to the squared ratio between the symmetry breaking scale and the Planck mass M p \ = 
G~ x / 2 . In general, we would expect this to be very small (e.g. on the order of 10 -6 for GUT scale strings), however, 
we will typically also investigate the solutions for higher values of k to understand the general pattern of solutions. 
The components of the energy-momentum tensor (in units of Xirjf) are 
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T t * = y^£j + u , Tl = -e\ + e 2 - £3 - £4 + £5 + u 



i=i 



Tg = £1 - e 2 - £3 - £4 + £5 + u , T* = £1 +e 2 +£3 -£4 - £5 + u , (11) 



where 



_ h' 2 f 2 _ h 2 P 2 _ 1 P' : 

£l - — + — ' 62 ' 63 ~ ~2a~T? ' (12) 



4= 2iV^ ' 5= 2iP ' U= 4 {h _1) + T /(/ - 2 l) + Y kf ■ (13) 
The three independent Einstein equations then read 

R\ = ~n{2e 4 + e 3 -u) , (14) 
Rg = n(2e 2 + e 3 +u) , (15) 
R\ = k (2£ 5 - £ 3 + u) (16) 

subject to the constraint (which is not independent) 

R r r = /t(2£!+£ 3 + u) . (17) 

The Euler-Lagrange equations which result from the variation of the action with respect to the matter fields are 

(LKNh'Y , „, x , „, hP 2 

± = h (h 2 - 1) + 73/1/ 2 + —r- , (18) 

LKN K L y ' 

(LKNf)' k 2 f u; 2 f 2 2 2 

ZiTiV = -^2-+73fe/ + 72/(/ -9), (19) 

L fNKP' 



NK V £ 



= aP/i' . (20) 
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Equations (fT4)) - (|T6| and (fT8|) - (l20|) have to be solved subject to appropriate boundary conditions. At r — the 
requirement of regularity and the fact that we would like string-like solutions leads to 

h(0) = , /'(0) = , P(0) = n , L(0) = , L'(0) = 1 , iV(0) = 1 , N'(0) = , K(0) = 1 , K'(0) = . (21) 

Finiteness of energy requires that 

h(oo) = 1 , /(oo) = , P(oo) = . (22) 

In our numerical calculations, we have solved the coupled system (|14 p -(|16 |> . (fT8")) - (l2T)| subject to the boundary 
conditions (|2"Tj) . (|22[). In all our calculations, we have checked that (|17p is fulfilled. 



B. Energy per unit length, tension and current 

In [24l ] the energy per unit length U and tension T have been used as macroscopic parameters to describe supercon- 
ducting string solutions in the f7(l)i OC ai x t/(l)giobai model in flat space-time and to investigate the stability of these 
objects. The corresponding expressions in curved space-time are a straightforward generalization of that and read 

-it 



U = J J V-hTl drd6 = 2tt / / L T\ dr , T = - J V-h T* drdO = -2vr / L T* dr, (23) 

where h corresponds to the determinant of the induced metric on the (t, z)-plane. It is also possible to define the 
so-called Tolman mass of these solutions which corresponds to the gravitational active mass (see e.g. [HI), but we 
would like to compare our results to the flat space-time limit studied in (T3 | and hence define T and U as given above. 
Furthermore, the Noether current associated to the unbroken U(l) symmetry reads 

jm = i (f*V"£ - CV£*) , (24) 

which has non-vanishing, dimensionless components 

jt = u W ' JZ = k i- (25) 

This current is covariantly conserved V M J p = 0, which implies <9 M (y / —gJ tl ) = and since \/—Q is independent of 
(t, z) we have J M = 0. Finally, we can also define the charge number density (in analogy to [Tj]) which reads 



J = 2n ^\J l J t + J z J z \ L dr . (26) 



In [15| a macroscopic stability condition for superconducting strings has been suggested. This requires that the 
propagation speeds of transverse (T) and longitudinal (L) perturbations, respectively, should both be real for the 
string to be stable. The squared propagation speeds are given by 

9 T dT . „ N 

c --u ' c - = ~du (27) 

and should hence both be positive. Since our definitions of T and U are covariantly defined expressions and the 
original work of Carter was formulated in a fully covariant way flH [l6| , we expect these criteria to also hold in the 
gravitating case and use them in the following to decide about the stability of gravitating superconducting strings. 



C. Asymptotic behaviour of the metric functions 



Far away from the string core where the matter functions have reached their vacuum values, we would expect that 
the metric functions have the behaviour of those of a general cylindrically symmetric vacuum space-time given by the 
Kasner metric which has the following form 



la / N 2b / x 2c 

.2 2 
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where a, b and c are real coefficients subject to the Kasner conditions: 

a + b + c = , a 2 + 6 2 + (c+l) 2 = 1 , (29) 

while r a can be thought of as the typical radius of the string. 7 is an additional parameter that determines the deficit 
angle of the space-time A. Note that the first Kasner condition can be used to eliminate c in the second one. We 
then get a second-order polynomial in b, which has real solutions if and only if a £ [— -|, l] . Then b is given by: 



1 - a± \/-3a 2 + 2a + 1 , 1 + a ± yJ-Za 2 + 2a + 1 



cF (30) 



2 ' 2 

so that = b T — 1. 

In [24| it was argued that with the identification 

a = -b = 2G(U - T) + 0(G 2 ) , c = 0{G 2 ) , 7 = 1 — 4G(U + T) (31) 

the Kasner metric would describe the outside of a general cosmic string with energy per unit length U and tension T 
under the assumption that the string is infinite, straight and has negligible width. The deficit angle of the space-time 
is then given by A = 4nG(U + T). 

For "standard" cosmic strings we have U — T and the metric describes a conical space-time with deficit angle 
A = 8irGU [2]. In this case, the space-time is locally flat and geodesies are just straight lines. Since the space-time 
has a deficit angle, gravitational lensing appears, but planetary, i.e. bound orbits are not possible. Note that this 
changes when taking the finite core of the cosmic string into account [2l[ . 

III. NUMERICAL RESULTS 

The solutions to the equations ( fP fll- flTrH) and (|18l) - (|2"0|) are only known numerically. We have solved these equations 
using the ODE solver COLSYS 1291. The solutions have relative errors on the order of 1CT 6 - 1(T 10 . In the following, 
we have restricted our analysis to the case n = 1 unless otherwise stated. 

A. General behaviour and Kasner coefficients 

In FigQ]we show the matter and metric functions of a solution for w = —0.3, 72 = 3, 73 = 2, q — 0.7, a — 0.1 and 
three different values of k. As is apparent, the matter functions vary only very little (the different cases are barely 
distinguishable on the plot), however, the metric functions change quite strongly. When increasing k the metric 
functions N(r) and K(r) deviate more strongly from their flat space-time values N(r) = K(r) = 1, while the slope of 
L(r) at large r decreases with increasing k signaling -as expected - an increase in the deficit angle of the space-time. 

The outside space-time of a superconducting string should be given by the Kasner space-time with a specific relation 
between the Kasner coefficients and U and T (see (|3T|) ) [24| . We observe that our numerical solutions indeed possess 
an asymptotic behaviour governed by the Kasner metric (|28p . Matching the space-time of the solution given in FigfT] 
at large r with the Kasner metric we find a = 0.00040092, b = -0.00040076, c = -0.00000054 for « = 0.001 and 
a = 0.04098961, b = -0.03937305, c = -0.00165716 for k = 0.1, respectively. 

In Figf2]we give the values of the Kasner coefficients in dependence on k(U — T) for a superconducting string with 
timclike current (the results are qualitatively similar in the spacelike case, this is why we do not present them here) 
and 72 = 3, 73 = 1.74, q = 0.69, a = 0.1. For small values of k we recover the behaviour given in (|3Tj) . For larger 
values of k the parameter c starts to deviate from zero, while a and b no longer depend linearly on k(U — T). We 
hence conclude that for small values of k(U — T) the description of the metric outside the superconducting string by 
a Kasner metric with the identification (|31l) is a very good approximation. However, with our techniques, we are also 
able to determine the behaviour of the metric functions inside the string core up to the string axis at r = 0. This is 
only possible when solving the equations of motion numerically. In the following, we want to discuss the macroscopic 
stability of gravitating superconducting strings. For this, we will need to integrate the matter and metric functions 
from r — to r = 00. It is hence crucial to know the matter and metric functions on the full interval re [0 : 00 [. 

B. Macroscopic stability 

In flat space-time it was found [l4| that in the spacelike regime w > the energy per unit length U and the tension 
T are always positive and hence c\ > 0. However, while the energy per unit length is an increasing function of w in 
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log(r+1) log(r+1) 



(a) matter functions P(r), h(r) : f(r) (b) metric function N(r) 




log(r+1) log(r+1) 



(c) metric function L(r) (d) metric function K(r) 



FIG. 1: The matter and metric functions for a gravitating superconducting string with w — —0.3, 72 — 3, 73 — 2, q — 0.7, a — 0.1 and different 

values of k. 



the spacelike regime, the tension T decreases only for w close to the chiral limit w = 0. For sufficiently small w it was 
hence found that the strings are stable with > 0, while for larger values cj* < 0. In the timelike regime, the energy 
per unit length and tension diverge at the approach of the phase frequency threshold which corresponds to the value 
of — w equal to the mass of the scalar boson. In the following, we have fixed 72 = 3, 73 = 2, a = 0.1, q = 0.7 unless 
otherwise stated. This choice of parameters fulfills all the requirements such that the local U(l) symmetry is broken 
and the global U(l) symmetry remains unbroken in the flat space-time limit. In the following, we will be interested 
in the way that the energy per unit length U, the tension T as well as the charge number density J change with k. 



In Fig|3]we plot U and T as function of w for w close to w = 0. For k = we recover the results given in 14 1. 

Note that we plot the dependence on w and not on sgn(w) \/]w|. We do not present our results in dependence on 
this latter quantity here, but have convinced ourselves that the plots look qualitatively similar to those presented in 
[l4j . We observe that the main features are still present for k 7^ 0. The value of w = w° r at which T = U in the 
spacelike regime corresponds to the value of w where /(0) = and hence f(r) = 0. We find that w° r increases with 
k, as do the corresponding values of U and T. This leads also to the observation that the range of w in the spacelike 
regime for which > increases since the minimal value of T appears at larger values of w. Hence, the coupling 
to gravity enhances the interval of w > in which strings are stable. This is also seen in Fig|4l where we plot T as 
function of U . Obviously, dT/dU < for w close to the chiral limit w = 0. We also give the charge number density 
J in FigEl This shows that this quantity vanishes at w° r as well as at w = 0. In the spacelike regime the charge 
number density first rises reaching a maximum at w max and then decreases again to zero at w[? r . We find that w max 
increases with k. In the timelike regime the charge number density rises strongly and the bigger k the bigger J for a 
given value of w. 

In the timelike regime we find in analogy with [T3 | that there exists a phase frequency threshold at which U — > +00 
and T — > — 00 for n = 0. This corresponds to the value of w at which scalar bosons can be produced. Interestingly, 
we observe that this phase frequency threshold seems to be absent when k ^ 0. This is shown in Fig|6l where we 
give U as function of w in the timelike regime. Clearly, for large values of — w the energy per unit length tends to a 
finite value for k^0. Interestingly, we find that the coupling to gravity can hence stabilize the strings. This is shown 
in FigO where we present U and T in the timelike regime for n = 0.3. In an interval close to w = we find that 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 

k-(U-T) 



FIG. 2: We show the value of the Kasncr coefficients a, b and c for gravitating superconducting string solutions with timelikc current and 
72 — 3, 73 — 1.74, q — 0.69, a — 0.1 in dependence on k(U — T). We give the coefficients for two different values of w. 




FIG. 3: We show the energy per unit length U (black) as well as the tension T (red) as function of w — k 2 — lj 2 for superconducting strings for 

different values of k and 72 — 3, 73 — 2, q — 0.7 and a — 0.1. 

both U and T are positive with U increasing and T decreasing such that cf^ > and > 0. Hence, the strings are 
stable. Decreasing w further leads to T becoming negative such that while c~[ > we now have c\ < and the strings 
are unstable. However, in contrast to the flat space-time case, where c\ < up to the phase frequency threshold, 
we observe that for sufficiently small w the tension becomes positive again. As such c\ > 0. Furthermore, for these 
values of w we find that U decreases, while T increases such that > and the strings are macroscopically stable. 
We conclude that there are hence two stable regions in the timelike regime for k^0. Hence, gravity can stabilize the 
strings for large values of the current. To get an idea how the parameter ranges in which the string becomes stable 
depend on k, we show the values of w where T = in dependence on n in FigJS] In between the two curves the 
tension T is negative and hence is negative. We observe that the value of w close to w = at which T vanishes 
and dT/dw < increases with increasing k. On the other hand, the value of w at which T vanishes and dT/dw > 
decreases with increasing k. Hence, the range w for which c\ < decreases with increasing k. Furthermore, our 
numerical results indicate that c\ > for all w if k is sufficiently large. For 72 = 3, 73 = 2, a = 0.1, q = 0.7 we find 
that this happens at k ~ 0.16. 
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FIG. 4: We show the tension T as function of the energy per unit length U for superconducting strings for different values of k and 72 — 3, 

73 = 2, q = 0.7 and a = 0.1. 




FIG. 5: We show the charge number density J as function of w — k — ui for superconducting strings for different values of k and 72 — 3, 

73 = 2, q = 0.7 and a = 0.1. 



IV. MOTION OF TEST PARTICLES 



In order to probe space-times it is crucial to understand how test particles move in these. We here consider 
structure-less, point-like particles that move on geodesies in the gravitational field of a superconducting string. If we 
consider the numerically given space-time in terms of the metric functions N, L and K we have to solve the geodesic 
equation numerically. This approach has been taken in [2l[ for standard cosmic strings and it was found that the 
finite core width of the string allows for additional bound orbits of massive test particles close to the string core. Since 
the string core and hence also the radius of the orbits are on the same order as the inverse scalar boson mass these 
cannot account for planetary orbits. On the other hand, the string could be "dressed" by massive particles and this 
could be important for gravitational wave emission. However, our numerical results indicate that the power of this 
emission would be rather small. This is why we do not investigate this further here. 

As shown above, we find that the metric outside the string core is well matched by the Kasner space-time. Since 
this space-time can be given analytically we are able to make some general statements about the geodesic motion. 

Since the observation of gravitational lensing by a cosmic string has been suggested as a possibility to detect these 
objects directly, we discuss the motion of massless test particles in the Kasner space-time in the following. We also 
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FIG. 6: We show U as function of — w for different values of k and 72 — 3, 73 — 2, q — 0.7 and a — 0.1. 




comment on the motion of massive test particles. 



A. The geodesic equation 



The Kasner metric has three obvious Killing vectors: dt, dg and d z . The three associated constants of motion are: 



2a / N 2b , / N 2c 

at „ I r \ dz 



where A is an affine parameter that can be identified with proper time for time-like geodesies. P t corresponds to the 
energy, P z to the linear momentum in the ^-direction and Pg to the angular momentum of a test particle moving in 
this space-time, respectively. The geodesic Lagrangian then reads 

1 dx^ dx v _ 1 
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T>0 



0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 

K 

FIG. 8: We show the region in the w-At domain for which T < (and hence < 0). Note that for k — the tension T > in the spacclike 

regime. This changes when considering the superconducting string in curved space-time for sufficiently small, but non-vanishing k. This 

is for 72 — 3, 73 — 2, q — 0.7 and a — 0.1. 

where e takes the value and 1 for massless and massive particles, respectively. Using the constants of motion (j3"2")l 
the equation of motion for r reads 



_ dt n dz _ de [dry 

Pt dX +P *dX +Pe dX + U) 



(34) 



which can be rewritten as follows 



2,i 



e I r 



4a 



P z \ 2 ( r x ! " 
~P t 



1 I r \ I Pi 



^r" \ r, 



= 



(35) 



The dynamics is thus similar to that of a classical point particle of unit mass and zero total energy subject to the 
effective potential: 



Kff(r) 



2a 

r \ e 



if V r n 



P,\ I r 



Pi 



i / 7-> \ 2 / \ 4a-2c 
1 I Pft\ I r 



jr 2 \P t 



(36) 



The equations for 9 and z are directly obtained from the constants of motion and read 



dz 
It 



Pt \r*. 



d9 
~dt 



J_Pi (r_ 

ryj.2 p t 



2(a-c) 



(37) 



B. The effective potential 



Some properties of the orbits can be deduced from the form of the effective potential. First note that V e g < in 
order for (|35|) to have solutions. 

In the following we will discuss the possible orbits in this space-time. We will define bound, i.e. planetary orbits 
as orbits on which particles move from a maximal finite radius r max to a minimal radius r m ; n > and back again in 
the x-y-pl&ne. Note that these orbits can extend to infinity in the z-direction for P z ^ 0. In contrast to that we have 
unbound orbits on which particles approach the string from infinity, move to a minimal radius r m i n and then return 
to infinity in the a;-y-plane. In both cases, it is also of interest to know whether r m i n can become zero, i.e. whether 
the particle can reach the string axis. In the following, we will call orbits that end at r m j n = "terminating orbits" . 
The possible orbits are summarized in Table I. 

In Figin] and FigfTUlwe show the effective potential V e s(r) for massive and massless particles, respectively, with 
P t = 2, P z = Pg = 0.2 in the space-time of a Kasner metric with r a = 1, 7 = 0.8 and different choices of a. The 
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n ffp nf r 




A 




terminating escape orbit 


B 1 




terminating orbit 


C 2 




bound orbit 


D 1 




escape orbit 



TABLE I: Types of possible orbits for tost particles moving in the static Kasncr space-time. 
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(a) 



effective potential V & ff(r) for b , c 



(b) 



effective potential V e ff(r) for b , c 



FIG. 9: We show the effective potential V e ff (r) for a massive test particle (e — 1) with Pt — 2, P z — Pq — 0.2 in a Kasncr space-time with 

r a — 1, 7 — 0.8 and different choices of a. Note that the Kasner conditions then fix the possible values of b and c. We show the effective 

potential for 6^, (left) as well as for b~ , c~ (right). 



choice of a fixes the value of b and c (see ([50)) ). We give the potential for both choices of b and c and note that the 
difference between the two choices is marginal. For general choices of a our results are summarized in Table II. These 
are: 

• If a < 0, there are no bound orbits since the potential has no local maximum and is negative at infinity. This 
is clearly seen in FigJH] and Fig llOl for a = —0.2. In addition, the minimal radius r m ; n is finite, except for a 
massless particle with vanishing P z and Pg. In this case, V c r (r) cx r 2a for r — > 0. Since a < 0, J£ oc ^ oc r~ a 
is integrable around r = and the particle thus reaches the string after a finite coordinate time. 

• For a = 0, there are no bound orbits. This corresponds either to the case of a standard cosmic string in which 
case the space-time is locally flat with b = b~ = c~ = c = or to b = b + = —c + = — c = 1 which is not physical. 
The minimal radius r mm is generally finite, except if P z = (for 6 + , c + ) or Pg = (for 6~, c~), in which case 
the potential is just a constant. 

• For a > and e — 1 (massive particle): the orbits are always bound if a < 2/3 and P z ^ (for b + , c + ) or 
a < 2/3 and Pg ^ (for b~ , c~). This is clearly seen in Fig|H]for a = 0.2 and a = 0.5, where the potential in 
both cases has positive values close to r = and tends to positive values at r — > oo. For other choices of a the 
orbits are bound terminating, i.e. end at the string axis at r = since the potential tends to zero from below 
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> and < 2/3 


> 2/3 
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type 


A 


D 


D 


D 


A 


D 


A 


B 


D 


C 


B 


c 


A 


B 


B 


B 



TABLE II: Types of orbits in the static Kasner space-time in dependence on the different parameters for the (b^ , c~^) solution. The table for 

(b~ , c~ ) would be the same with P z and Pe exchanged. 
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at r — > (see the potential for a = 0.8 in Figj9j). When the minimal radius vanishes, the dominant term near 
r = is at least r 2a . So, 4^ = (D{r~ a ). Since a < 1, this is integrable around r = and the time needed to 
reach the string is finite (except if a = 1). 

• For a > and e = (massless particle): the behaviour around r — is exactly the same as in the massive case 
(since the removed term is negligible), but the orbits are unbound for a < 2/3 when choosing Pq — for (b + , 
c + ) or P z = for (&~, c~). This is seen in Fig llOl This means that for appropriate choices of the parameters, 
we find that bound orbits for massless test particles exist. 





(a) 



effective potential V e ff(r") for c+ 



(b) 



effective potential V^ff(r) for b , c 



FIG. 10: We show the effective potential V c ff(r) for a massless test particle (e — 0) with Pt — 2. P z — P$ — 0.2 in a Kasncr space-time with 

r CT — 1, 7 — 0.8 and different choices of a. Note that the Kasner conditions then fix the possible values of b and c. We show the effective 

potential for (&~*~, c~*~) (left) as well as for (b~~ , c~) (right). 



We remark that exchanging Pq and P z does not change the qualitative form of the potential, provided we also 
exchange c + ) with c~). In fact, the equation of motion for r is invariant under 

(6+,c+)o(6- O , P z ^-^= (38) 

so that we can restrict our attention to one of the two cases. 

Studying the asymptotic properties of V e s is a priori not enough to rule out bound orbits, or to say that there is 
only one bound orbit (the potential could cross zero several times in the intermediate region.) However, we find that 
nothing new arises from a more detailed study. 

Note that the case of interest in the context of cosmic strings with additional structure corresponds to a > 0, a as 
0, b = b~ « —a, c = c~ w 0. 



C. Numerical results 



The equations of motion are solved numerically using a second-order symplectic integrator, which avoids numerical 
dissipation effects. We studied the domain of existence and the dependence of the maximal and minimal radius on 
the Kasner coefficients and constants of motion, as well as the light deflection using Maple and Mathematica. 



1. Domain of existence 

Given the parameters a, b and c in the Kasner metric, a natural question is what values of 7, P<, P z and Pe allow 
for solutions of the equations of motion. A necessary and sufficient condition is that the potential must be negative 

somewhere. For a < or a > 2/3, the dominant term at zero or infinity respectively is — \ (f-J • Therefore, 
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solutions always exist. For a £ ]0, 2/3 [, the domain in which solutions exist can be computed numerically. For 
massless particles, since r a can be taken equal to one by rescaling r and 7 can be taken to one by rescaling Pg, the 
function r i-» V e g(r) depends only on two parameters, for instance on P z /Pt and Pg/Pt for r a — 1, 7 = 1. Moreover, 
because of the aforementioned duality, we can restrict our attention to the (b~,c~) solution. For massive particles, 
Pt becomes a physically relevant parameter. The domain in which solutions exist tend towards the one obtained for 
e = in the limit Pt — > 00 and shrink when Pt decreases. For a = 0, we found analytically that solutions exist if 
and only if P^ + e < P? for the (6 + ,c + ) case or — + e < P t 2 for the (b~,c~) case (in units r a = 1). For a — 2/3, 








(a) P t = 1 (b) Pt = 2 

FIG. 11: Wc show the domain of existence of solutions to the geodesic equation in the Pg-P s -planc in a general Kasncr space-time with v a — 1, 
7—1 and different values of a for massless test particles with Pt — 1 (left) and Pt — 2 (right). Note that solutions to the geodesic 

equation exist below and to the left of the curves. 



2. Examples of orbits 

The most important difference to the case T = U in which geodesies are just straight lines because the space-time is 
locally flat is that we can have bound orbits for massive and massless test particles in the general Kasner space-time. 
In FigfT5]we show bound orbits for a massive and massless test particle, respectively. The bound orbit of a massive 
test particle shows the typical perihelion shift of orbits in curved space-time. The bound orbit of the massless test 
particle also possesses this perihelion shift, however looks qualitatively rather like a spiral. We observe that this in 
related to the choice of a, which is quite big. For a small the orbit of a massless particle looks qualitatively similar to 
that of a massive one. 

In Fig[T3] we show escape orbits for a massive and massless test particle, respectively. In both cases, the test 
particles get deflected by the string, which is related to the presence of the deficit angle of the space-time. We observe 
that the massive test particle experiences a stronger deflection for the same values of all parameters as the massless 
one. 



3. Minimal and maximal radius of bound orbits 

If the momenta are held fixed, the minimal and maximal radii of the bound orbits are decreasing functions of a. 
The maximal radius generally goes to infinity when a — > + , while the minimum radius remains finite. This was 
to be expected since the orbits are unbound with finite minimal radius for a < (let alone the very special case 
Pg = P z — 0). For a > 2/3 and fixed Pg, Pt, we find terminating bound orbits. The maximal radius decreases with 
P z , going to zero at infinity. It remains finite as P z — > 0, which was also expected from the above analysis. For 
< a < 2/3, there are two critical values of P z above which the equations of motion can not be satisfied for (b + , c + ) 
or (b~ , c~). The minimal and maximal radius merge at this point, and move apart from each other when we decrease 
P z . The merging of the minimal and maximal radius corresponds to a local minimum of the effective potential. If we 
consider the (& + ,c + ) case, the maximal radius remains finite when P z — > 0, while the minimum one goes to zero. For 
the (b~ , c~) case, it is r m i n which remains finite while r max goes to infinity if e = 0. For a < 0, we find no maximal 
radius, and the minimal one increases with P z . In Fig|14l we show the minimal and maximal radius of the bound 
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(a) massive test particle, Pt — 2, P z — 0, Pq — 0.5 , (b) massless test particle, P± — 2.0, P z — 0.1, P$ — 0.03 , 
a = 0.3, b = -0.2266281297, c = -0.07337187027 and a = 0.8, b = 0.512311, c = -1.31231 and 7 = 0.8 

7 = 0.5 



FIG. 12: We show a bound orbit for a massive (left) and a masslcss (right) test particle in a Kasncr spacc-timc. 




(a) massive test particle (b) massless test particle 

FIG. 13: We show escape orbits for a massive test particle (left) and a masslcss test particle (right) for P± — 2. P z — 0.1. Pq — 0.2 , a — 0, 

b — b— and c — c_ and 7 — 0.5 in a Kasncr spacc-timc. 

orbits for a Kasner space-time with a = 0.3, b = -0.2266281297, c = -0.073371870270 and 7 = 0.5 in dependence on 
Pg. The test particle has Pt = 2 and three different values of P z . The r m ; n and r max curves meet at the circular orbit. 
We find that the interval of Pg for which bound orbits exist increases when decreasing P z for a fixed P t . Moreover, 
we observe that the radius of the circular orbit decreases with increasing P z . In addition, the circular orbit appears 
at smaller angular momentum Pg when increasing P z . 
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0.5 1 p 1.5 2 2.5 

8 

FIG. 14: The value of the minimal radius r m i n and of the maximal radius r max of bound orbits of massive test particles (e — 1) with Pt — 2 in a 
Kasncr space-time with a — 0.3, b~ , c~ , 7 — 0.5 arc shown in dependence of Pq for three different values of P z . 



4- Light deflection 

For an unbound orbit, we can define the deviation angle as: 

AS = 2 r " dr -„=^r 1(1.) ~" -A= ir - „ (39) 
Jr min dr 7 Pt J rmin r 2 \rj v ; 

where r m i n is the minimal distance from the string to the particle. 

If all the other parameters are fixed, A6 is a decreasing function of 7. It goes to infinity when 7 — >■ and is equal 
to ^-^= — 1^ 7r in conical space-time. (This is true in a coordinate system in which the metric has the Kasncr form 

and 9 goes from to 2ir. A geometrical argument shows that it reduces to the usual 2tt(1 — ^Pf) in a locally flat 
coordinate system.) 

For a > 0, the (b~,c~) case gives an unbound orbit for massless particles with vanishing P z . We find that the 
deviation angle increases with a. 

The deviation angle also seems to behave logarithmically in Pg in the limit Pg — > 00. For Pg — > we have A8 —> —n. 
This is related to the repulsive effect of the string: if the angular momentum vanishes, the particle just goes back to 
where it comes from. 



V. SUMMARY AND DISCUSSIONS 



In this paper we have studied superconducting string solutions of the E/(l)i OC ai x ^(l)global model in curved space- 
time. For small values of the ratio between the symmetry breaking scale and the Planck mass we find that the metric 
outside the string core matches the Kasner metric very well and the Kasner coefficients can be given in terms of the 
energy per unit length U and the tension T as originally suggested in 24]. In order to be able to decide about the 
macroscopic stability of these objects using the stability conditions given by Carter [H [ij] we have to know the 
metric functions on the full interval r £ [0 : 00 [. We have hence integrated the full coupled system of differential 
equations numerically and computed the energy per unit length and tension. We find that the coupling to gravity can 
stabilize the strings with large charge number density for sufficiently large values of the ratio between the symmetry 
breaking scale and the Planck mass and that in general the phase frequency threshold is absent in curved space-time. 
As such, the energy per unit length and tension never diverge, but tend to more or less constant values for larger 
values of the charge number density. 

We have also studied the motion of test particles in the general Kasner space-time and find that the fact that the 
energy per unit length is non-equal to the tension can lead to bound orbits of massive and massless test particle. 
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The question is then whether these orbits are of interest for astrophysical or cosmological applications. We find that 
the radius of these orbits is too large to be of interest and that e.g. the gravitational wave emission from a particle 
moving on a bound orbit around a cosmic string would be far too small to be detectable. This - in turn - means 
that the assumption of an infinitely thin cosmic string that is often used in simulations of string networks is a valid 
assumption and that the additional structure on the string does not have a big influence on the results. 
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